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Abstract. We examine the properties of the atmo- 
spheres of extrasolar giant planets at orbital distances 
smaller than 0.1 AU from their stars. We show that these 
"51 Peg b- like" planets are rapidly synchronized by tidal 
interactions, but that small departures from synchronous 
rotation can occur because of fluid-dynamical torques 
within these planets. Previous radiative-transfer and evo- 
lution models of such planets assume a homogeneous at- 
mosphere. Nevertheless, we show using simple arguments 
that, at the photosphere, the day-night temperature dif- 
ference and characteristic wind speeds may reach ~ 500 K 
and ~ 2 km sec -1 , respectively. Substantial departures 
from chemical equilibrium are expected. The cloud cov- 
erage depends sensitively on the dynamics; clouds could 
exist predominantly either on the dayside or nightside, 
depending on the circulation regime. Radiative-transfer 
models that assume homogeneous conditions are there- 
fore inadequate in describing the atmospheric properties 
of 51Pegb-like planets. We present preliminary three- 
dimensional, nonlinear simulations of the atmospheric cir- 
culation of HD209458b that indicate plausible patterns for 
the circulation and generally agree with our simpler esti- 
mates. Furthermore, we show that kinetic energy produc- 
tion in the atmosphere can lead to the deposition of sub- 
stantial energy in the interior, with crucial consequences 
for the evolution of these planets. Future measurements 
of reflected and thermally-emitted radiation from these 
planets will help test our ideas. 

Key words: extrasolar planets, 51 Peg, HD209458, tides, 
atmospheric dynamics, interior structure, giant planets 



1. Introduction 

The discovery of extrasolar planets has led to a growing 
list of work devoted to modeling their atmospheres (Bur- 
rows et al. 1997; Seager and Sasselov 1998, 2000] ; Gouken- 



leuque et al. 2000 



the similar brown dwarfs, for which theoretical models 
now reproduce the observations well, even in the case of 
low-tempe ratur e objects (T c ff ~ 100 K or less) (e.g. Mar- 
ley et al. |l996t Al lard et al. |l997| ; Lieber t et al. |2000 ; 



Geballe et al. |2001| ; Schweitzer et al. |2001| to cite only a 
few). However, an important feature of extrasolar plan- 
ets is their proximity to a star: The irradiation that they 
endure can make their atmospheres significantly different 
than those of isolated brown dwarfs with the same ef- 
fective temperature. This has been shown to profoundly 
alter the atmos pheric vertical temperatur e pro file (Sea- 
ger a nd Sas selov 1998 ; Goukenleuque et al. 200C ; Barman 
et al. 2001). We will show that most importantly, it also 



2001). While no spec- 



Barman et al. 

tra of these objects have yet been measured, one might 
be encouraged by the successes obtained in the case of 



affects the horizontal temperature distribution and atmo- 
spheric chemistry so that the models calculated thus far 
may fail to provide an adequate description of the atmo- 
spheres of the most intensely irradiated planets. Advec- 
tion has never previously been considered, but it can play 
a major role for the composition, temperature, spectral 
appearance and evolution of extrasolar planets. 

We will focus on the extrasolar planets for which irra- 
diation is the most important: 51 Peg b- like planets, which 
we henceforth dub "Pegasi planets" and define as gas gi- 
ants orbiting solar-type stars at less than 0.1 AU. Their 
importance is demonstrated by the fact that they orbit 
nearly 1% of stars surveyed so far and constitute 27% of 
currently-known extrasolar giant planets. They are also 
more easily characterized by the transit method than are 
other planets, as indicated by the discovery of the transit- 
ing gas giant HD209458b (Charbonneau et al. [2000 ; Henry 
et al. 200C). In the preceeding paper (Guillot and Show- 
man 2001, hereafter Paper I), we showed how the atmo- 
spheric boundary condition governs the evolution of Pe- 
gasi planets. We also advocated that an additional source 
of energy is needed to explain the radius HD209458b, and 
that this would most likely be provided by the down- 
ward transport and subsequent dissipation of kinetic en- 
ergy with a flux of ~ 1% of the absorbed stellar energy. 
In this paper, we use the temperature profiles obtained 
in Paper I to evaluate the dynamical state of Pegasi- 
planet atmospheres, and discuss how dynamics may in- 
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flucnce the cloud abundance, chemical composition, and 
thermal state (all of which will be amenable to observation 
in the near future). We then present preliminary numer- 
ical simulations of the circulation that indicate plausible 
circulation patterns and show how downward propagation 
of kinetic energy from the atmosphere to the interior can 
occur. 

After reviewing the expected interior structure (Sec- 
tion 2), we begin in Section 3 with the problem of tides: Pe- 
gasi planets have been predicted to rotate synchronously 
(Guillot et al. 199E), implying that they always present 
the same face toward the star. We argue however that 
dynamical torques may maintain the interior in a non- 
synchronous rotation state, which has important impli- 
cations for understanding atmospheric processes. In sec- 
tion 4, we discuss the probable wind speeds, day-night 
temperature differences, and flow geometries, including 
both order-of-magnitude arguments and our numerical 
simulations. A summary of the results is provided in sec- 
tion 5. 

2. Interior structure 

Fully self-consistent models including both the atmo- 
spheres and interiors of Pegasi planets do not yet exist. 
To obtain first-cut estimates of the expected tempera- 
ture profiles, we therefore use the same approach as in 
Paper I: we choose two extreme evolution models ("hot" 
and "cold") that match the radius of HD209458b (but at 
different ages and hence intrinsic luminosities). These evo- 
lution calculations only pertain to pressures larger than 3 
to 10 bars (depending on the boundary condition), and so 
we extend the profiles to lower pressure using radiative- 
transfer calculations for intensely-irradiated planets from 
Marley (personal communication) and Barman (2001). 
The calculations are not strictly appropriate to the case of 
HD209458b, but provide us with reasonable estimates of 
the expected structure (and as we will see our qualitative 
conclusions are not sensitive to the type of profile chosen). 
The resulting "hot" and "cold" models that we will use 
hereafter are depicted in Figure [l]. Interestingly, Paper I 
shows that these two models are representative of a rel- 
atively wide variety of models of HD209458b, including 
those with energy dissipation. 

It should be noted that the high temperatures can trig- 
ger the dissociation of the hydrogen molecule, which is in- 
dicated in Fig. |l] by a local maximum of the heat capacity 
c p . This effect is important for our purposes because it 
implies that stellar heat can be stored in hot regions and 
reclaimed in colder regions by molecular recombination. 
Other molecules (e.g. H 2 0) can undergo dissociation, but 
this will be neglected due to their small abundances. 

Evolution models of Pegasi planets show that two in- 
terior regions exist: a radiative zone (including the atmo- 
sphere) that extends down to pressures of 100 to 800 bar, 
and a deeper convective core. The intrinsic luminosity of 
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Fig. 1. Left: Temperature profiles for the "hot" model 
(thin line) and "cold" model (thick grey line) at times 
(5.37 and 0.18 Ga, respectively) when the models match 
HD209458b's measured radius. The diamonds indicate 
the radiative/convective boundary. The terms "hot" and 
"cold" are chosen to indicate the models' relative tem- 
peratures in the region from ~ 1-300 bars, which is the 
important region of the atmosphere for affecting the evo- 
lution. See Paper I for details. The discontinuity in the 
profiles' slopes is unphysical; it results from the fact that 
atmospheric models assume a much higher intrinsic lumi- 
nosity than predicted by evolution models. Right: Dimen- 
sionless mean specific heat per particle in the atmosphere. 
The local maxima in c p are due to the dissociation of the 
H2 molecule. 



the planet is about 10,000 times smaller than the energy 
absorbed from the star. If energy dissipation in the outer 
layers is large (see Paper I), another convective zone can 
appear at the levels where dissipation is the largest. This 
however requires a very high dissipation of ~ 10% of the 
incoming heat flux. We therefore chose to only consider 
the simple radiative/convective scenario, as depicted in 
Fig. |. 

We split the planet into an "atmosphere" dominated 
by stellar heating, with possible horizontal temperature 
inhomogeneities, and an "interior" including the con- 
vective core, for which inhomogeneities should be much 
smaller. The crux of the problem is to understand how 
the heat absorbed on the day side and near the equator is 
redistributed by winds and/or rotation to the night side 
and to the poles. 

3. Synchronization of Pegasi planets 

It has been shown that the tides raised by the star on 
Pegasi planets should ra pidly drive them int o sync hronous 
rotati on (G uillot et al. |l996[ Marcy et al. |l997t Lubow 
et al. 1997). This can be seen by considering the time 
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Fig. 2. Conjectured dynamical structure of Pegasi plan- 
ets: At pressures larger than 100-800 bar, the intrinsic 
heat flux must be transported by convection. The con- 
vective core is at or near synchronous rotation with the 
star and has small latitudinal and longitudinal tempera- 
ture variations. At lower pressures a radiative envelope is 
present. The top part of the atmosphere is penetrated by 
the stellar light on the day side. The spatial variation in 
insolation should drive winds that transport heat from the 
day side to the night side (see text). 



scale to tidally despin the planet (Goldreich , 
Hubbard g|8§): 



Soter 1966 



<syn 
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(1) 



where Q, R, M, a, ui and ui s are the planet's tidal dissi- 
pation factor, radius, mass, orbital semi-major axis, rota- 
tional angular velocity, and synchronous (or orbital) an- 
gular velocity. M* is the star's mass, and G is the gravi- 
tational constant. Factors of order unity have been omit- 
ted. A numerical estimate for HD209458b (with cu equal 
to the current Jovian rotation rate) yields a spindown 
time T syn ~ 3Q years. Any reasonable dissipation factor 



Q (see Marcy et al. 1997 Lubow et al. 1997| ) shows that 
HD209458b should be led to synchronous rotation in less 
than a few million years, i.e., on a time scale much shorter 
than the evolution timescale. Like other Pegasi planets, 
HD209458b is therefore expected to be in synchronous ro- 
tation with its 3.5-day orbital period. 

Nevertheless, stellar heating drives the atmosphere 
away from synchronous rotation, raising the possibility 
that the interior's rotation state is not fully synchronous. 
Here, we discuss (1) the energies associated with the 
planet's initial transient spindown, and (2) the possible 
equilibrium states that could exist at present. 



3.1. Spindown energies 

Angular momentum conservation requires that as the 
planet spins down, the orbit expands. The energy dissi- 
pated during the spindown process is the difference be- 
tween the loss in spin kinetic energy and the gain in orbital 
energy: 



E = -4- ( \k 2 MR 2 uj 2 - Va 2 ^ 
dt \ 2 2 



(2) 



where k is the dimensionless radius of gyration (k 2 = 
I /MR 2 , I being the planet's moment of inertia). The 
time derivative is negative, so E, the energy dissipated, 
is positive. The orbital energy is the sum of the planet's 
gravitational potential energy and orbital kinetic energy 
and is negative by convention. The conservation of angu- 
lar momentum implies that the rate of change of lj s is 
constrained by that on u>: 



— (Ma 2 u s + k 2 MR 2 uj) = 0. 
dt ' 



(3) 



The fact that the planetary radius changes with time may 
slightly affect the quantitative results. However, since r syn 
is so short, it can be safely neglected in this first-order 
estimate. R being held constant, it is straightforward to 
show, using Kepler's third law, that: 



E 



-k 2 MR 2 (uj -lo s )lu. 



(4) 



(Note that ill is negative, and so E is positive.) 

The total energy dissipated is E « k 2 MR 2 {uj s — w) 2 /2, 
neglecting variation of the orbital distance. Using the mo- 
ment of inertia and rotation rate of Jupiter (fc 2 — 0.26 
and lu = 1.74 x 10~ 4 sec" 1 ), we obtain for HD209458b 
£«4x 10 41 erg. If this energy were dissipated evenly 
throughout the planet, it would imply a global tempera- 
ture increase of 1400 K. 

By definition of the synchronization timescale, the dis- 
sipation rate can be written: 



E = 



k 2 MR 2 (u -lu s ) 



(5) 



' syn 



With Q of 10 5 , a value commonly used for Jupiter, 
3 x 10 5 years and the dissipation rate is then 



' syn 



10 erg sec , or 35,000 times Jupiter's intrinsic lumi- 



nosity. Lubow et al. (1997) have suggested that dissipation 
in the radiative zone could exceed this value by up to two 
orders of magnitude, but this would last for only ^100 
years. 

The thermal pulse associated with the initial spin- 
down is large enough that, if the energy is dissipated in 
the planet's interior, it may affect the planet's radius. It 
has previously been argued (Burrows et al. 2000] ) that 
Pegasi planets must have migrated inward during their 
first 10 7 years of evolution; otherwise, they would have 
contracted too much to explain the observed radius of 
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HD209458b. But the thermal pulse associated with spin- 
down was not included in the calculation, and this extra 
energy source may extend the time over which migration 
was possible. 

Nevertheless, it seems difficult to invoke tidal synchro- 
nization as the missing heat source necessary to explain 
HD209458b's present (large) radius. High dissipation rates 
are possible if T syn is small, but in the absence of a mech- 
anism to prevent synchronization, E would drop as soon 
as t > T syn . The most efficient way of slowing the planet's 



'sy: 

contraction is then to invoke r, 



syn 



10 



m . 



In that 
l 



years 

case, the energy dissipated becomes E ~ 10 24 erg sec 
which is two orders of magnitude smaller than that neces- 
sary to significantly affect the planet's evolution (Paper I; 
Bodenheimer et al. 2001). For the present-day dissipation 
to be significant, an initial rotation rate 10 times that of 
modern-day Jupiter would be needed. But the centripetal 
acceleration due to rotation exceeds the gravitational ac- 
celeration at the planet's surface for rotation rates only 
twice that of modern-day Jupiter, so this possibility is 
ruled out. Furthermore, such long spindown times would 
require a tidal Q of ~ 10 9 -10 10 , which is ~ 10 4 times the Q 
values inferred for Jupiter, Uranus, and Neptune from con 
straints on th eir satellites' orbits (Peale 
Murray 1992, Tittemore and Wisdom 



1999 , Banfield and 
1989). Dissipation 



of the energy due to transient loss of the planet's initial 
spin energy therefore cannot provide the energy needed to 
explain the radius of HD209458b. 

Another possible source of energy is through circular- 
ization of the orbit. Bodenheimer et al. (2001) show that 



the resulting energy dissipation could reach 10 26 ergsec -1 
if the planet's tidal Q is 10 6 and if a hypothetical com- 
panion planet pumps HD209458b's eccentricity to values 
near its current observational upper limit of 0.04. If such a 
companion is absent, however, the orbital circularization 
time is ~ 10 8 years, so this source of heating would be neg- 
ligible at present. Longer circularization times of 10 9 -10 10 
years would allow the heating to occur until the present- 
day, but its magnitude is then reduced to 10 25 erg sec -1 
or lower, which is an order of magnitude smaller than the 
dissipation required. 



3.2. The equilibrium state 

The existence of atmospheric winds implies that the atmo- 
sphere is not synchronously rotating. Because dynamics 
can transport angular momentum vertically and horizon- 
tally (including the possibility of downward transport into 
the interior), the interior may evolve to an equilibrium ro- 
tation state that is asynchronous. Here we examine the 
possibilities. 

Let us split the planet into an "atmosphere" , a part of 
small mass for which thermal effects are significant, and an 
"interior" encompassing most of the mass which has min- 
imal horizontal thermal contrasts. Suppose (since r syn is 
short) that the system has reached steady state. Consider 



two cases, depending on the physical mechanisms that de- 
termine the gravitational torque on the atmosphere. 

The first possibility is that the gravitational torque 
on the atmosphere pushes the atmosphere away from 
synchronous rotation (i.e., it increases the magnitude of 
the atmosphere's angular momentum measured in the 
synchronously-rotating reference frame) as has been hy- 



pothesized for Venus (Ingersoll and Dobrovolskis 1978 
Gold and Soter 1969Q . To balance the torque on the atmo- 



sphere, the interior must have a net angular momentum 
of the same sign as the atmosphere (so that both either 
super- or subrotate). In Fig. [| superrotation then corre- 
sponds to a clockwise flux of angular momentum, whereas 
subrotation corresponds to the anticlockwise scenario. 

The second possibility is that the gravitational torque 
on the atmosphere tends to synchronize the atmosphere. 
This possibility may be relevant because, on a gas-giant 
planet, it is unclear that the high-temperature and high- 
pressure regions would be 90° out of phase, as is expected 
on a terrestrial planet. If the interior responds sufficiently 
to atmospheric perturbations to keep the deep isobars 
independent of surface meteorology (an assumption that 
seems to work well in modeling Jupiter's cloud-layer dy- 
namics), then high-pressure and high-temperature regions 
will be in phase, which would sweep the high-mass re- 
gions downwind and lead to a torque that synchronizes the 
atmosphere. Furthermore, if a resonance occurs between 
the tidal frequency and the atmosphere's wave-oscillation 
frequency, the resulting gravitational torques also act to 
synchronize the atmosphere (Lubow et al. 1997). In this 
case, the interior and atmosphere have net angular mo- 
menta of opposite signs. Depending on the sign of the 
atmosphere/interior momentum flux, this situation would 
correspond in Fig. || to either the clockwise or anticlock- 
wise circulation of angular momentum cases. 

In both cases above, the gravitational torque on the at- 
mosphere arises because of spatial density variations asso- 
ciated with surface meteorology, which are probably con- 
fined to pressures less than few hundred bars. Therefore, 
this torque acts on only a small fraction of the planet's 
mass. In contrast, the gravitational torque on the gravi- 
tational tidal bulge affects a much larger fraction of the 
planet's mass. It is not necessarily true that the torque on 
the atmosphere is negligible in comparison with the torque 
on the interior, however, because the gravitational tidal 
bulge is only slightly out of phase with the line-of-sight 
to the star (e.g., the angle is 10~ 5 radians for Jupiter). In 
contrast, the angle between the meteorologically-induced 
density variations and the line of sight to the star could be 
up to a radian. Detailed calculations of torque magnitudes 
would be poorly constrained, however, and will be left for 
the future. 

Several mechanisms may act to transfer angular mo- 
mentum between the atmosphere and interior. Kelvin- 
Helmholtz shear instabilities, if present, would smooth 
interior-atmosphere differential rotation. On the other 
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Fig. 3. Angular momentum flow between orbit, interior, 
and atmosphere for a Pegasi planet in steady state. Ar- 
rows indicate flow of prograde angular momentum (i.e., 
that with the same sign as the orbital angular momen- 
tum) for two cases: Anticlockwise: Gravitational torque 
on atmosphere is retrograde (i.e., adds westward angular 
momentum to atmosphere). For torque balance, the gravi- 
tational torque on the interior must be prograde (i.e., east- 
ward). These gravitational torques must be balanced by 
fluid-dynamical torques that transport retrograde angular 
momentum from atmosphere to interior. Clockwise: Grav- 
itational torque on atmosphere is prograde, implying a 
retrograde torque on the interior and downward transport 
of prograde angular momentum from atmosphere to inte- 
rior. Atmosphere will superrotate if gravitational torques 
push atmosphere away from synchronous (as on Venus). 
It will subrotate if gravitational torques synchronize the 
atmosphere (e.g., gravity- wave resonance; cf. Lubow et al. 
(1997)). 



hand, waves transport angular momentum and, because 
they act nonlocally, often induce differential rotation 
rather than removing it. The quasi-biennial oscillation on 
Earth is a classic example (see, e.g., Andrews et al. 1987, 
Chapter 8). Atmospheric tides (large-scale waves forced 
by the solar heating) are another example of such a wave; 
on Venus, tides play a key role in increasing the rotation 
rate of the cloud-level winds. Finally, vertical advection 
may cause angular momentum exchange between the at- 
mosphere and the interior of Pegasi planets, because the 
angular momentum of air in updrafts and downdrafts need 
not be the same. 

A simple estimate illustrates the extent of nonsyn- 
chronous rotation possible in the interior. Suppose that 
the globally-averaged flux of absorbed starlight is F 9 , 
which is of order 10 s ergs -1 cm -2 for Pegasi planets near 
0.05 AU, and that the globally-averaged flux of kinetic 
energy transported from the atmosphere to the interior is 
T)F a , where i] is small and dimensionless. If this kinetic 
energy flux is balanced by dissipation in the interior with 
a spindown timescale of T syn , then the deviation of the 
rotation frequency from synchronous is 



3 syn 



k 2 M 



1/2 



This estimate is an upper limit for the globally averaged 
asynchronous rotation rate because it assumes that kinetic 
energy transported into the interior has angular momen- 
tum of a single sign. The global-average asynchronous ro- 
tation could be even lower if angular momenta transported 
downward in different regions have opposite signs. 

The spindown time is uncertain and depends on the 
planet's Q according to Eq. (1). To date, only one es timat e 
for the Q of a Pegasi planet exists (Lubow et al. 1997 ), 
which suggests Q ~ 100, at least during the early stages 
of spindown. Orbital constraints on the natural satellites 
of the giant planets imply that, at periods of a few days, 
the tidal Q of Uranus, Neptune, and Jupiter arc of order 
~ 1 5 (T ittemo re an d Wisdom |l989j , Banfield and Mur- 
ray 1992 , Peale 1999| ), which suggests a spindown time of 
a few xlO 5 years (Eq. 1). Although the mechanisms that 
determine these planets' Q values remain uncertain, likely 
possibilities include friction in the solid inner core (Der- 
mott 1979) or overturning of tidally- forced waves in the 



(6) 



fluid envelope (Ioannou and Lindzen 1993, Houben et al. 
[2001 ). Both mechanisms are possible for Pegasi planets, 
and the wave-dissipative mechanism should be more ef- 
fective for Pegasi planets than for Jupiter because of the 
thicker stable (radiative) layer in the former (Houben et 
al. 2001). These calculations suggest that spindown times 
of ~ 10 5 ~10 6 years are likely, although larger values can- 
not be ruled out. 

As discussed in Section 4, experience with planets in 
our solar system suggests that atmospheric kinetic energy 
is generated at a flux of 10 _2 F (> , and if all of this energy 
enters the interior, then r\ ~ 10 -2 . Using a spindown time 
of 3 x 10 5 years implies that u> — oj s ~ 2 x 10 -5 scc -1 , 
which is comparable to the synchronous rotation fre- 
quency. The implied winds in the interior are then of or- 
der ~ 2000 msec -1 . Even if ?y is only 10 -4 , the interior's 
winds would be 200 msec -1 . The implication is that the 
interior's spin could be asynchronous by up to a factor 
of two, depending on the efficiency of energy and momen- 
tum transport into the interior. If spindown times of ~ 10 6 
years turn out to be serious under- or overestimates, then 
greater or lesser asynchronous rotation would occur, re- 
spectively. 



4. Atmospheric Circulation: Possible regimes and 
influence on evolution 

4-1. Basic Parameter Regime 

Rotation plays a central role in the atmospheric dynamics 
of Pegasi planets, and HD209458b in particular. The ratio 
of nonlinear advective accelerations to Coriolis accelera- 
tions in the horizontal momentum equation is it/ -1 L -1 
(called the Rossby number), where u is the mean hori- 
zontal wind speed, / = 2Vt sin <f) is the "Coriolis param- 
eter" (e.g. Holton |l992| , pp. 39-40), fl is the rotational 
angular velocity, <j) is latitude, and L is a characteristic 
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length scale. Rossby numbers of 0.03-0.3 are expected for 
winds of planetary-scale and speeds ranging from 100- 
1000 msec -1 . In Section 3 we showed that modest asyn- 
chronous rotation may occur, in which case the Rossby 
number could differ from this estimate by a factor up 
to ~ 2. These estimates suggest that nonlinear advective 
terms are small compared to the Coriolis accelerations, 
which must then balance with the pressure-gradient accel- 
erations. (The Rossby number could be ~ 1 if the winds 
reach 2000-3000 msec -1 , which we show below is proba- 
bly the maximum allowable wind speed.) 

The zonality of the flow can be characterized by the 
Rhines' wavenumber, kp ~ (/3/u) 1 / 2 , where (3 is the 
derivative of / with northward distance (Rhines 1975). 
The half- wavelength implied by this wavenumber, called 
the Rhines' scale Lp, provides a reasonable estimate for 
the jet widths on all four outer planets in our solar system 
(Cho and Polvani 1996; see Table 1). For HD209458b, the 
Rhines' scale is ~ 1.5 x 10 10 (w/1000 msec -1 ) 1 / 2 cm, which 
exceeds the planetary radius if the wind speed exceeds 
about 400 msec -1 . 

Another measure of horizontal structure is the Rossby 
deformation radius (Gill 1982, p. 205), L D ~ NH/f, 
where H is the scale height and N is the Brunt- Vaisala 
frequency (i.e., the oscillation frequency for a vertically 



displaced air parcel; Holton 1992, p. 54). At pressures of a 
few bars, the temperature profiles calculated for irradiated 
extrasolar giant planets by Goukenlouque et al. (2000) 
suggest N ~ 0.0015 sec -1 . With a scale height of 700 km, 
the resulting deformation radius is 40,000 km. In contrast, 
the deformation radii near the tropopause of the Jupiter, 
Saturn, Uranus, and Neptune are of order 2000 km (Table 
!)■ 

The estimated Rhines' scale (for winds of > 
500 msec -1 , which we show later are plausible speeds) 
and deformation radius of Pegasi planets are similar to the 
planetary radius, and they are a larger fraction of the plan- 
etary radius than is the case for Jupiter, Saturn, Uranus, 
and Neptune (Table 1). This fact suggests that eddies may 
grow to hemispheric scale in the atmospheres of Pegasi 
planets and that, compared with the giant planets in our 
solar system, the general circulation hot Jupiters may be 
more global in character. Unless the winds are extremely 
weak, Pegasi planets are unlikely to have > 10 jets as do 
Jupiter and Saturn. 

An upper limit on the atmospheric wind speed can 
be derived from shear-instability considerations. We as- 
sume that no zonal winds are present (u(Pq) = 0) in 
the convective core, a consequence of synchronization by 
tidal friction. The build-up of winds at higher altitudes in 
the radiative envelope is suppressed by Kelvin-Helmholtz 
instabilities if the shear becomes too large. This occurs 



when the Richardson number becomes smaller than 1/4 
(cf. Chandrasekhar 1961), i.e. when 



Ri = 



N 2 1 

(du/dzf < 4' 



(7) 



In the perfect gas approximation, assuming a uniform 
composition, 



N 



H 



(V ad - V T ). 



(8) 



where 1Z is the universal gas constant divided by the mean 
molecular mass, c p is the specific heat, H = IZT/g is the 
pressure scale height, V a d = lZ/c p is the adiabatic gradi- 
ent, Vt = d In T/dlnp, T is temperature, and p is pres- 
sure. The hypothesis of uniform composition is adequate 
despite hydrogen dissociation because the timescales in- 
volved are generally much longer than the dissociation 
timescales. It can be noted that a shear instability could 
appear at smaller Richardson numbers in the presence 
of efficient radiative diffusion (Zahn 1992; Maeder 1995). 
This possibility will not be examined here. 

The maximal wind speed at which Kelvin-Helmholtz 
instabilities occur can then be derived by integration of 
Eq. <§: 



c(P) 



[KT(V ad -V T )} 1/2 d\np, 



(9) 



Pi I 



The value of w max thus derived at P ~ 1 bar is of the or- 
der of 3000 m s -1 , to be compared to the winds of Jupiter, 
Saturn, Uranus and N eptune which reach 100-500 ms -1 
(Ingersoll ct al. 1995). Our estimate for Pegasi planets, 



assuming the winds are measured in the synchronously- 
rotating frame, may be uncertain by a factor of ~ 2 due 
to the possibility of a non-synchronously-rotating inte- 
rior. In comparison, the expected speed of sound at the 
tropopause of HD209458b is ~ 2400 msec -1 . 

A characteristic timescale for zonal winds to redis- 
tribute temperature variations over scales similar to the 
planetary radius R then stems from r zona i > -R/u max . 

The radiative heating timescale can be estimated by 
a ratio between the thermal energy within a given layer 
and the layer's net radiated flux. In the absence of dy- 
namics, absorbed solar fluxes balance the radiated flux, 
but dynamics perturbs the temperature profile away from 
radiative equilibrium. Suppose the radiative equilibrium 
temperature at a particular location is T ra d and the actual 
temperature is T ra d + AT. The net flux radiated towards 
outer space is then 4aT^ d AT and the radiative timescale 
is 



P 



^rad 



g 4crT 3 ' 



(10) 



where a is the Stefan-Boltzmann constant. 

Figure [| shows estimates of T zona i and r ra d for 
HD209458b calculated using the temperature profiles from 
the "hot" (thin line) and "cold" (thick grey line) models 
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Table 1. Rhines' and deformation lengths for giant planets 





u 


R 


n 




L D 


Jet width 




(msec - ) 


(10 7 m) 


(lO^sec- 1 ) 


(10 7 m) 


(10 7 m) 


(10 7 m) 


Jupiter 


50 


7.1 


1.74 


1.0 


0.2 


~ 1 


Saturn 


200 


6.0 


1.6 


1.9 


0.2 


~ 2 


Uranus 


300 


2.6 


1.0 


2.0 


0.2 


~ 2 


Netpune 


300 


2.5 


1.1 


1.9 


0.2 


~ 2 


HD209458b 


? 


10 


~ 0.2 


15(u/lkmsec _1 ) 1/2 


~ 4 


? 



Note. Ld calculated at the tropopause using N = 0.01 sec for Saturn, Uranus, and Neptune, and N = 0.02 sec for Jupiter. 



from Section 2. The zonal timescale is estimated by calcu- 
lating the maximum wind speed that can exist as a func- 
tion of pressure given the static stability associated with 
each model, while the radiative time is calculated using 
the temperatures and heat capacities shown in Fig. [I]. At 
pressures exceeding 0.1 bar, radiation is slower than the 
maximal advection by zonal winds, but by less than one 
order of magnitude. The consequent day/night tempera- 
ture difference ATday-night to be expected is: 



ATriav- 



day — night 



1 - e" 



a< 1 



(11) 



where AT ra d is the day-night difference in radiative equi- 
librium temperatures. Rough estimates from Fig. [| sug- 
gest that r zona i/r r ad ~ 0.3 at 1 bar, implying that 
AT day _ night /AT rad -0.3. If AT rad = 1000 K, this would 
imply day-night temperature differences of 300 K at 1 bar. 
Values of ATday-night even closer to AT ra d are likely given 
the fact that slower winds will lead to an even more ef- 
fective cooling on the night side and heating on the day 
side. 

The small radiative time scale implies that, for the 
day-night temperature difference to be negligible near the 
planet's photosphere, atmospheric winds would have to 
be larger than the maximum winds for the onset of shear 
instabilities. 




4 5 6 

Log(Timescale) [s] 



7 -4 -3 -2 -1 
Log(q/c p ) [K/s] 



Fig. 4. Left: Characteristic time scales as a function of 
pressure level. r zona i is the minimal horizontal advection 
time (dashed). r ra d is the timescale necessary to cool/heat 
a layer of pressure P and temperature T by radiation alone 
assuming the gas is optically thick (solid). (The plotted 
radiative timescale is an underestimate at pressures less 
than about 0.3 bars because the atmosphere becomes op- 
tically thin at those pressures.) For each case, the thin 
black and thick grey lines correspond to the hot and cold 
models from Section 2 and Figj^. Right: Approximate cool- 
ing/heating rate as a function of pressure (see Eq. 111). 



4-2. Nature of the circulation 

An understanding of the horizontal temperature difference 
and mean wind speed is desirable. It is furthermore of in- 
terest to clarify the possible geometries the flow may take. 
For Pegasi planets, we envision that the dominant forcing 
is the large-scale day- night heating contrast, with minimal 
role for moist convection. This situation differs from that 
of Jupiter, where differential escape of the intrinsic heat 
flux offsets the solar heating contrast and moist convection 
plays a key role. For Pegasi planets, such an offset between 
intrinsic and stellar fluxes cannot occur, because the in- 
trinsic flux is 10 4 times less than the total flux. The fact 
that the day-night heating contrast occurs at hemispheric 
scale — and that the relevant dynamical length scales for 
Pegasi planets are also hemispheric (Section 4.1) — in- 
creases our confidence that simple analysis, focusing on 
the hemispheric-scale circulation, can provide insight. 



Because the Rossby number is small, the dominant 
balance in the horizontal momentum equation at mid- 
latitudes is between the Coriolis force and the pressure- 
gradient force (geostrophic balance). We adopt the prim- 
itive equations, which are the standard set of large-scale 
dynamical equations in a stably-stratified planetary at- 
mosphere. When differentiated with pressure (which is 
used here as a vertical coordinate), this balance leads to 



the well-known thermal wind equation (e.g. Holton 1992 
75): 



P 



<9v 
d\np 



-TZk x V H T 



(12) 



where v is the horizontal wind, p is pressure, Vh is the 
horizontal gradient, T is temperature, and k is the unit 
upward vector. Assuming the interior winds are small 



8 



A. P. Showman and T. Guillot: Dynamics of Pegasi-planet Atmospheres 



compared to the atmospheric winds, this implies that, to 
order-of-magnitudc, 



TZ 



z Alnp 



(13) 



where |v| and AThoriz are the characteristic wind speed 
and horizontal temperature difference in the atmosphere 
and Amp is the difference in log-pressures from bottom 
to top of the atmosphere. Non-sychronous rotation of the 
interior would alter the relation by changing the value of 
/, but the uncertainty from this source is probably a factor 
of two or less (see Eq. (6) and subsequent discussion). 

The thermodynamic energy equation is, using pressure 
as a vertical coordinate (e.g. Holton 1992, p. 60), 



dT 



— h v • VhT — w- 



H 2 N 2 
TZp 



(14) 



where t is time, u> = dp/dt is vertical velocity, q is the 
specific heating rate (ergg -1 sec -1 ), and we have assumed 
an ideal gas. 

A priori, it is unclear whether the radiative heating 
and cooling is dominantly balanced by horizontal advec- 
tion (second term on left of Eq. (14)) or vertical advec- 
tion (third term on the left). To illustrate the possibilities, 
we consider two endpoint scenarios corresponding to the 
two limits. In the first scenario, the radiation is balanced 
purely by horizontal advection: zonal winds transport heat 
from dayside to nightside, and meridional winds transport 
heat from equator to pole. In the second scenario, the radi- 
ation is balanced by vertical advection (ascent on dayside, 
descent on nightside). 

Consider the first scenario, where horizontal advection 
dominates. Generally, we expect v and VhT to point in 
different directions, so to order of magnitude their dot 
product equals the product of their magnitudes. (If the 
direction of VhT is independent of height and no deep 
barotropic flow exists, then at mid-latitudes one could ar- 
gue that winds and horizontal pressure gradients are per- 
pendicular to order Ro. However, the existence of either 
asynchronous rotation or variations in the orientations of 
VhT with height will imply that winds and VhT are not 
perpendicular even if Ro << 1. Because asynchronous ro- 
tation and variation in the directions of VhT with height 
are likely, and because, in any case, v and VhT will not be 
perpendicular near the equator, we write the dot product 
as the product of the magnitudes.) An order-of-magnitude 
form of the energy equation is then 



|v|AT horiz _ q_ 



R 

The solutions are 

1/2 



q TZA In p 



(15) 



(16) 



AXho 



R 



f 



c p TZA In p 



1/2 



(17) 



A rough estimate of the heating rate, q/c p , results from 
the analysis in Section 4.1: 



q _ AaT 3 ATg 

Cp PCp 



(18) 



where AT is the characteristic magnitude of the differ- 
ence between the actual and radiative equilibrium tem- 
peratures in the atmosphere. The heating rate depends on 
the dynamics through AT. We simply evaluate the heat- 
ing rate using Eq. ( |l§| ) with AT « T/2. The results are 
shown in Fig. ||. 

The key difficulty in applying the equations to Pe- 
gasi planets is the fact that q/c p depends on pressure, 
and AThoriz probably should too, but Eqs. (16) and (17) 
were derived assuming that AThoriz is constant. We can 
still obtain rough estimates by inserting values of q/c p at 
several pressures. At ~ 50-100 bars, where q/c p ~ 10 -4 - 
10 -3 K sec -1 , we obtain temperature differences and wind 
speeds of 50-150 K and 200-600 msec -1 . At 1 bar, where 
q/cp reaches 10 -2 K sec -1 , the estimated temperature con- 
trast and wind speed is ~ 500 K and ~ 2000 msec -1 . 
The estimates all assume / w 3 x 10 -5 sec -1 , TZ = 
3500 J kg K -1 , and Alnp«3. 

Later we show that the equations successfully pre- 
dict the mean wind speeds and temperature differences 
obtained in numerical simulations of the circulation of 
HD209458b. This gives us confidence in the results. 

Now consider the second scenario, where vertical ad- 
vection (third term on the left of Eq. (14)) balances the ra- 
diative heating and cooling. The magnitude of u) can be es- 
timated from the continuity equation. Purely geostrophic 
flow has zero horizontal divergence, so uj is roughly 



W ~ R °~k~ P 



(19) 



Strictly speaking, this is an upper limit, because the non- 
linear terms comprising the numerator of the Rossby num- 
ber contain some terms (e.g., the centripetal accelera- 
tion) that are divergence-free. Substituting this expres- 
sion into the energy equation, using Eq. (13), and setting 
Ro = |v|//i? implies that 



nhW 



1/2 



fR 2 (q_l 
NHA\np \c p TZ 



1/2 



(20) 
(21) 



Numerical estimates using H = 700 km, N ~ 
0.0015 sec -1 , and q/c p = 10 -4 -10 -3 K sec -1 (appropriate 
to 50-100 bars) yield temperature differences and speeds 
of 80-250 K and 300-900 msec -1 , respectively. Here we 
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have used the same values for R and / as before. A heating 
rate of 10 _2 Ksec -1 , appropriate at the 1 bar level, yields 
|v| - 3000msec- 1 and AT horiz ~ 800 K. These values are 
similar to those obtained when we balanced heating solely 
against horizontal advection. 

We can estimate the vertical velocity under the as- 
sumption that vertical advection balances the heating. Ex- 
pressing the vertical velocity w as the time derivative of 
an air parcel's altitude, and using w = —w/pg, where p is 
density, yields 



Using q/c p = 10 -2 Ksec -1 and H = 700 km implies that 
w ~ 20 msec -1 near 1 bar. If this motion comprises the 
vertical branch of an overturning circulation that extends 
vertically over a scale height and horizontally over ~ 10 10 
cm (a planetary radius), the implied horizontal speed re- 
quired to satisfy continuity is ~ 3000 msec" 1 , consistent 
with earlier estimates. 

The numerical estimates, while rough, suggest that 
winds could approach the upper limit of ~ 3000 msec -1 
implied by the shear instability criterion. This compari- 
son suggests that Kclvin-Helmholtz shear instabilities may 
play an important role in the dynamics. 

The likelihood of strong temperature contrasts can 
also be seen from energetic considerations. The differ- 
ential stellar heating produces available potential energy 
(i.e., potential energy that can be converted to kinetic 
energy through rearrangement of the fluid; Peixoto and 
Oort 1992, pp. 365-370), and in steady-state, this po- 
tential energy must be converted to kinetic energy at 
the rate it is produced. This requires pressure gradi- 
ents, which only exist in the presence of lateral ther- 
mal gradients. A crude estimate suggests that the rate 
of change of the difference in gravitational potential <f> 
between the dayside and nightside on isobars caused by 
the heating is ~ 1Z(q/c p )A Inp. The rate per mass at 
which potential energy is converted to kinetic energy by 
pressure-gradient work is V'Vij$, which is approximately 
|v|7?.AThorizAlnp/i?. Equating the two expressions sug- 
gests, crudely, that AThoriz ~ (g/ c p )( J R/|v|), which implies 
AThoriz ~ 500 K using the values of q/c p and |v| near 1 
bar discussed earlier. 

One could wonder whether a possible flow geometry 
consists of dayside upwelling, nightside downwclling, and 
simple acceleration of the flow from dayside to nightside 
elsewhere (leading to a flow symmetrical about the sub- 
solar point). The low-Rossby- number considerations de- 
scribed above argue against this scenario. A major com- 
ponent of the flow must be perpendicular to horizontal 
pressure gradients. Direct acceleration of wind from day- 
side to nightside is still possible near the equator, how- 
ever. This would essentially be an equatorially-confined 
Walker- type circulation that, as described earlier, would 



lie cquatorward of the geostrophic flows that exist at more 
poleward latitudes. 

4-3. Numerical Simulations of the Circulation 

To better constrain the nature of the circulation, we 
performed preliminary three-dimensional, fully-nonlinear 
numerical simulations of the atmospheric circulation of 
HD209458b. For the calculations, we used the Explicit 
Planetary Isentropic Coordinate, or EPIC, model (Dowl- 
ing et al. 1998). The model solves the primitive equations 
in spherical geometry using finite-difference methods and 
isentropic vertical coordinates. The equations are valid in 
stably-stratified atmospheres, and we solved the equations 
within the radiative layer from 0.01 to 100 bars assuming 
the planet's interior is in synchronous rotation with the 
3.5-day orbital period. The radius, surface gravity, and 
rotation rate of HD209458b were used (10 8 m, 10 msec -2 , 
and f2 = 2.1 x 10 -5 scc -1 , respectively). 

The intense insolation was parameterized with a sim- 
ple Newtonian heating scheme, which relaxes the temper- 
ature toward an assumed radiative-equilibrium tempera- 
ture profile. The chosen radiative-equilibrium temperature 
profile was hottest at the substellar point (0° latitude, 0° 
longitude) and decreased toward the nightside. At the sub- 
stellar point, the profile's height-dependence was isother- 
mal at 550 K at pressures less than 0.03 bars and had 
constant Brunt- Vaisala frequency of 0.003 sec -1 at pres- 
sures exceeding 0.03 bars, implying that the temperature 
increased with depth. The nightside radiative-equilibrium 
temperature profile was equal to the substellar profile mi- 
nus 100 K; it is this day-night difference that drives all the 
dynamics in the simulation. The nightside profile was con- 
stant across the nightside, and the dayside profile varied 
as 

Tdaysido = Tnightsido + AT rad COS OL (23) 

where Tdaysido is the dayside radiative-equilibrium tem- 
perature at a given latitude and longitude, T n i g htsido is 
the nightside radiative-equilibrium temperature, AT ra d = 
100 K, and a. is the angle between local vertical and the 
line-of-sight to the star. For simplicity, the timescale over 
which the temperature relaxes to the radiative-equilibrium 
temperature was assumed constant with depth with a 
value of 3 x 10 5 sec. This is equal to the expected radiative 
timescale at a pressure of about 5 bars (Fig. 4, left). 

The Newtonian heating scheme described above is, 
of course, a simplification. The day-night difference in 
radiative-equilibrium temperature, 100 K, is smaller than 
the expected value. Furthermore, the radiative timescale 
is too long in the upper troposphere (~ 0.1-1 bar) and too 
short at deeper pressures of ~ 10-100 bars. Nevertheless, 
the net column-integrated heating per area produced by 
the scheme (Wm -2 ) is similar to that expected to occur 
in the deep troposphere of HD209458b and other Pegasi 



10 



A. P. Showman and T. Guillot: Dynamics of Pegasi-planet Atmospheres 



planets, and the simulation provides insight into the cir- 
culation patterns that can be expected. _ 

The temperature used as the initial condition was of 
isothermal at 500 K at pressures less than 0.03 bars and^ 
had a constant Brunt- Vaisala frequency of 0.003 sec -1 at ^ 
pressures exceeding 0.03 bars. There were no initial winds. ^ 
The simulations were performed with a horizontal resolu- 2> 
tion of 64 x 32 with 10 layers evenly spaced in log-pressure. 

Figure 5-8 show the results of such a simulation. De- 
spite the motionless initial condition, winds rapidly de- 
velop in response to the day-night heating contrast, reach- 
ing an approximate steady state after ~ 400 Earth days."™ 
Snapshots at 42 and 466 days are shown in Figs. 5 and^ 
6, respectively. In these figures, each panel shows pressure ■S 
on an isentrope (greyscale) and winds (vectors) for three 
of the model layers corresponding to mean pressures of S 
roughly 0.4, 6, and 100 bars (top to bottom, respectively). 
The greyscale is such that, on an isobar, light regions are 
hot and dark regions are cold. 

The simulation exhibits several interesting features. 
First, peak winds exceed 1km sec -1 , but despite these"™ 
winds, a horizontal temperature contrast is maintained. 2_ 
The dayside (longitudes —90° to 90°) is on average hotter % 
than the nightside, but dynamics distorts the temperature 3 
pattern in a complicated manner. Second, an equatorial jet « 
develops that contains most of the kinetic energy. The jet 
initially exhibits both eastward and westward branches 
(Fig. 5) but eventually becomes only eastward (Fig. 6), 
and extends from —30° to 30° in latitude. Third, away 
from the equator, winds develop that tend to skirt paral- 
lel to the temperature contours. This is an indication that 
geostrophic balance holds. (Because the thermal structure 
is independent of height to zcroeth order and we have as- 
sumed no deep barotropic flow, horizontal pressure and 
temperature gradients are parallel.) Nevertheless, winds 
are able to cross isotherms near the equator, and this is 
important in setting the day-night temperature contrast. 
As expected from the order-of-magnitude arguments in 
Section 4.1, the jets and gyres that exist are broad in 
scale, with a characteristic width of the planetary radius. 

In the simulation, the day-night temperature differ- 
ence (measured on isobars) is about 50 K. This value de- 
pends on the adopted heating rate, which was chosen to 
be appropriate to the region where the pressure is tens 
of bars. Simulations that accurately predict the day-night 
temperature difference at 1 bar will require a more realis- 
tic heating-rate scheme; we will present such simulations 
in a future paper. 

The temperature patterns in Figs. 5 and 6 show that 
Earth-based infrared measurements can shed light on the 
circulation of Pegasi planets. In Fig. 6, the superrotating 
equatorial jet blows the high-temperature region down- 
wind. The highest-temperature region is thus not at the 
substellar point but lies eastward by about 60° in longi- 
tude. The maximum and minimum temperatures would 
thus face Earth before the transit of the planet behind 
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-100 100 




-100 100 

Longitude [deg] 
Fig. 5. EPIC simulation of atmospheric circulation on 
HD209458b after 42 days. Panels depict pressure on an 
isentrope (greyscale) and winds (vectors) for three model 
levels with mean pressures near 0.4, 6, and 100 bars (top to 
bottom, respectively). From top to bottom, the maximum 
wind speeds are 937, 688, and 224 msec -1 , respectively, 
and the greyscales span (from dark to white) 0.31-0.49 
bars, 5.6-7.8 bars, and 92-113 bars, respectively. Substel- 
lar point is at 0° latitude, 0° longitude. 



and in front of the star, respectively. On the other hand, 
if a broad westward jet existed instead, the maximum 
and minimum temperatures would face Earth after the 
transits. Therefore, an infrared lightcurve of the planet 
throughout its orbital cycle would help determine the di- 
rection and strength of the atmospheric winds. 

In the simulation, the intense heating and cooling 
causes air to change entropy and leads to vertical mo- 
tion, as shown in Fig. 7. Light regions (Fig. 7) indicate 
heating, which causes ascent, and dark regions indicate 
cooling, which causes descent. Air is thus exchanged be- 
tween model layers. This exchange also occurs across the 
model's lowermost isentrope, which separates the radia- 
tive layer from the convective interior (Fig. 7, bottom). 
Because upgoing and downgoing air generally have differ- 
ent kinetic energies and momenta, energy and momentum 
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Fig. 6. Simulation results for HD209458b at 466 days (af- 
ter a steady state has been reached). As in Fig. 5, panels 
depict pressure on isentropes (greyscale) and winds (vec- 
tors) for three model levels with mean pressures near 0.4, 
6, and 100 bars (top to bottom, respectively). From top 
to bottom, the maximum wind speeds are 1541, 1223, and 
598 msec -1 , respectively, and the greyscales span (from 
dark to white) 0.32-0.51 bars, 5.6-8.1 bars, and 90-127 
bars, respectively. Substellar point is at 0° latitude, 0° 
longitude. 



-100 100 

Longitude [deg] 
Fig. 7. Additional simulation results at 466 days. Ar- 
rows are identical to those in Fig. 6, but here greyscale 
is vertical velocity dO/dt, where 9 is potential tempera- 
ture. Light regions are heating (i.e., 9 is increasing) and 
dark regions are cooling (9 is decreasing). From top to 
bottom, the greyscales span (from dark to white) —0.003 
to 0.003 K sec" 1 , -0.002 to 0.002 K sec" 1 , and -0.001 to 
0.001 K sec -1 , respectively. The lowermost panel separates 
the convective interior from the radiative region, and the 
implication is that mass exchange can happen across this 
interface. 



can thus be exchanged between the atmosphere and inte- 
rior. 

Figure 8 indicates how the steady state is achieved. 
The build-up of the mass- weighted mean speed (top panel) 
involves two timescales. Over the first 10 days, the mean 
winds accelerate to 110 msec -1 , and the peak speeds ap- 
proach 1km sec -1 . This is the timescale for pressure gra- 
dients to accelerate the winds and force balances to be es- 
tablished. The flow then undergoes an additional, slower 
(~ 300 day) increase in mean speed to 300 msec -1 , with 
peak speeds of ~ 1.5kmsec _1 . This timescale is that 
for exchange of momentum with the interior, across the 
model's lowermost isentrope, to reach steady state, as 
shown in Fig. 8, middle. 



In the simulation, kinetic energy is transported from 
the atmosphere into the interior at a rate that reaches 
2500 Wm -2 (Fig. 8, bottom). This downward transport of 
kinetic energy is ^ 1% of the absorbed stellar flux and is 
great enough to affect the radius of HD209458b, as shown 
in Paper I. Although the simulation described here does 
not determine the kinetic energy's fate once it reaches the 
convective interior, we expect that tidal friction, Kelvin- 
Hclmholtz instabilities, or other processes could convert 
these winds to thermal energy. The exact kinetic energy 
flux will depend on whether a deep barotropic flow ex- 
ists. Furthermore, potential and thermal energy are also 
transported through the boundary, and pressure work is 
done across it. Our aim here is not to present detailed di- 
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agnostics of the energetics, but simply to point out that 
energy fluxes that are large enough to be important can 
occur. We are currently conducting more detailed simu- 
lations to determine the sensitivity of the simulations to 
a deep barotropic flow, and we will present the detailed 
energetics of these simulations in a future paper. 

The evolution of angular momentum (Fig. 8, middle) 
helps explain why the circulation changes from Fig. 5 to 6. 
After 42 days (Fig. 5), eastward midlatitude winds have 
already developed that balance the negative equator-to- 
pole temperature gradient. But at this time, the atmo- 
sphere's angular momentum is still nearly zero relative 
to the synchronously-rotating state, so angular momen- 
tum balance requires westward winds along some parts of 
the equator. (Interestingly, in these regions, the temper- 
ature increases with latitude, as expected from geostro- 
phy.) The circulation involves upwelling on the dayside 
and equatorial flow — both east and west — to the night- 
side, where downwelling occurs. This circulation is similar 
in some respects to the Walker circulation in Earth's at- 
mosphere. After 466 days, however, the atmosphere has 
gained enough angular momentum that the equatorial jet 
is fully superrotating. Some westward winds exist at high 
latitudes, but they are weak enough that the net angular 
momentum is still eastward. 

Application of the arguments from Section 4.2 pro- 
vides a consistency check. The mean heating rate in the 
simulations is about 1.7 x 10 -4 Ksec -1 . Inserting this 
value into Eqs. (16-17), we obtain ATh or iz = 70 K and 
|v| = 240 msec -1 . These compare well with the day- night 
temperature difference and mass-weighted mean speed of 
50 K and 300 msec -1 that are obtained in the simulation. 
(For the estimate, we used / = 3 x 10 -5 sec -1 , 1Z = 
3500 J kg -1 K -1 , R = 10 8 m, and Ap = 3.) If Eqs. (20- 
21) are adopted instead, using values of N — 0.003 sec -1 
and H = 300 km, we obtain temperature differences and 
speeds of 130 K and 470 msec -1 . Compared to the sim- 
ulation, these estimates are too high by a factor of ~ 2, 
but are still of the correct order of magnitude. An impor- 
tant point is that the maximum speed in the simulation 
substantially exceeds that from the simple estimates; the 
estimates are most relevant for the mean speed. 
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Fig. 8. Top: Mass-weighted mean wind speed within 
the model domain (~ 0.01-100 bar) during the sim- 
ulation of HD209458b shown in Figs. 5-7. Middle: 
Globally-averaged zonal angular momentum per area mi- 
nus that of the synchronously-rotating state. (That is, 
RA~ X J u cos (j) dm, where A is the planet's area, u is east- 
ward speed, (j) is latitude, and dm is an atmospheric mass 
element.) Bottom: Globally-averaged flux of kinetic en- 
ergy across the model's bottom isentrope (at ~ 100 bars), 
which is the interface between the radiative layer and the 
convective interior in this model. 



4-4- Vertical motion and clouds 

As shown in the appendix, hot and cold regions on the 
planet may have distinct chemical equilibrium composi- 
tions. Advection of air between these regions plays a key 
role for the (disequilibrium) chemistry in the atmosphere. 
Cloud formation critically depends on whether vertical 
motions dominate over horizontal motions. This in turns 
affects the albedo and depth to which stellar light is ab- 
sorbed. 

Clouds form when the temperature of moving air 
parcels decreases enough for the vapor pressure of con- 
densible vapors to become supersaturated. If dayside heat- 



ing and nightside cooling are balanced only by horizontal 
advection, air flowing from dayside to nightside remains 
at constant pressure but decreases in temperature (hence 
entropy) because of the radiative cooling. Nightside con- 
densation may occur, and particle fallout limits the total 
mass of the trace constituent (vapor plus condensates) in 
this air. As the air flows onto the dayside, its temperature 
increases and any existing clouds will sublimate. This sce- 
nario therefore implies that the dayside will be cloud-free. 

In an alternate scenario, dayside heating and nightside 
cooling are balanced by vertical advection, at least near 
the sub- and antisolar points. Because entropy increases 
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with altitude in a stably-stratified atmosphere, this sce- 
nario implies that ascent occurs on the dayside and de- 
scent on the nightside. Clouds may therefore form on the 
dayside, increasing the albedo and decreasing the pressure 
at which stellar light is absorbed. 

In reality, both horizontal and vertical advection are 
important, and the real issue is to determine the rela- 
tive contribution. We are pursuing more detailed numeri- 
cal simulations to address this issue. 



5. Conclusions 

We examined the structure and dynamics of the atmo- 
spheres of Pegasi planets and concluded that strong day- 
night temperature contrasts (300 K or more) are likely to 
occur near the level where radiation is emitted to space 
(~ 1 bar). These temperature variations should drive a 
rapid circulation with peak winds of 1km sec -1 or more. 
Force-balance arguments suggest that the mean midlati- 
tude winds are eastward, but the equatorial winds could 
blow either east or west. Depending on the dynamical 
regime, the cloud coverage and consequently the radia- 
tive absorption of the incoming stellar flux will be very 
different. 

Preliminary numerical simulations show that, while 
the dayside is generally hotter than the nightside, the spe- 
cific temperature distribution depends on the dynamics. 
In our simulations, a broad superrotating jet develops that 
sweeps the high-temperature regions downwind by about 
1 radian. This implies that the greatest infrared flux would 
reach Earth ~ 10-15 hours before the occultation of the 
planet behind the star. On the other hand, if a subro- 
tating jet existed, the greatest infrared flux would reach 
Earth after the occultation. Measurements of the infrared 
lightcurve of Pegasi planets will therefore constrain the 
direction and magnitude of the wind. 

The simulations also produce a downward flux of ki- 
netic energy across the ~ 100 bar surface equal to ~ 1% 
of the absorbed stellar flux. Although the simulation did 
not explicitly include the convective interior, we surmise 
that a substantial fraction of this kinetic energy would 
be converted to thermal energy by Kelvin-Helmholtz in- 
stabilities and other processes. As discussed in Paper I, 
"standard" evolution models explain HD209458b's radius 
only when the atmosphere is assumed to be unrealisti- 
cally hot, but addition of an internal energy source allows 
a more realistic model to match the observed radius. The 
downward transport and subsequent dissipation of kinetic 
energy described here is a promising candidate. Boden- 
heimer et al. (2001) suggested an alternative — that the 
internal heating could be provided by tidal circularization 
of an initially eccentric orbit. The difficulty, as Boden- 
heimer et al. were careful to point out, is that the tidal 
heating is a transient process in the absence of a detected 
close, massive companion capable of exciting the planet's 
eccentricity. In contrast, the mechanism proposed here can 



last throughout the star's lifetime. Nevertheless, to fully 
determine the feasibility of the mechanism, more detailed 
numerical simulations are required (in particular, to test 
the dependence of the energetics on the possible existence 
of winds in the interior) . We will present such simulations 
in a future paper. 

Upcoming observations are likely to provide key in- 
formation within the decade. Several spacecraft missions 
(either proposed or accepted) and dedicated ground pro- 
grams will observe extrasolar planets. Measurement of 
starlight reflected from these planets may allow the albedo 
to be estimated. Because the star-planet-Earth angle 
changes throughout the planet's orbit, crude informa- 
tion on the scattering properties of the atmosphere (e.g., 
isotropic versus forward scattering) may be obtainable. 
Asymmetries in the reflected flux as the planet approaches 
and recedes from the transit could give information on 
the differences of albedo near the leading and trailing ter- 
minators, which would help constrain the dynamics. Fi- 
nally, transit observations of Pegasi planets using high 
resolution spectroscopy should in the near future yield 
constraints on the atmospheric temperature, cloud/haze 
abundance, and winds (Brown 2001; see also Seager & 
Hubbard et al. 



Sasselov 2000 



2001). If these measure- 



ments are possible during the ingress and egress, i.e., the 
phases during which the planets enters and leaves the stel- 
lar limb, respectively, asymmetries of the planetary signal 
should be expected and would indicate zonal heat advec- 
tion at the terminator. The duration of these phases being 
limited to less than 10 minutes, it is not clear that this ef- 
fect is observable with current instruments. 
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Appendix A: Day-night temperature variations 
and chemical composition 

Here we explore the implications of day-night temperature 
variations on the abundances of a condensing species. Sup- 
pose that on the day side, a given condensate (e.g. iron) 
condenses at a pressure level P^ ay and a corresponding 
temperature Td a y If the temperature on the night side 
Tmght ls smaller, at what pressure -P* ight will condensa- 
tion take place? 

Assuming an ideal gas, we write the Clausius- 
Clapeyron equation as 



dlnp 
d\nT 



= /3 



(A.l) 



where p is the partial pressure of saturation of the con- 
densing species and f3 = L/1ZT is the ratio of the latent 
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heat of condensation L to the thermal energy TZT. For 
most condensing species of significance here, (3 » 10. We 
furthermore assume that j3 is independent of T and P, 
which introduces only a slight error in our estimates. 

By definition, on the day side, the saturation abun- 
dance of the condensing species, x — p/P is maximal and 
equal to x* at P = P£ ay . On the other hand, the night 
side temperature is lower and the abundance becomes: 



ln:r(P day ) = lux* - /31n(T day /r nigh t). 



(A.2) 



In order to reach condensation, i.e. x = x* . on e has to 
penetrate deeper into the atmosphere. Eq. ( A.l) implies 
that 



din: 



= /3V T - 1, 



dlnP 

and hence, on the night side, 

lnx(P) = lnx(P d * ay ) + 09V T - 1) hi(P/P d * ay ), 



(A.3) 



(A.4) 



assuming that Vj is constant. Using eq. (A.2), one obtains 
the condensation pressure on the night side: 



P* 

night 

p* 

May 



T 



day 



T, 



night 



/3/(/3V T -l) 



(A.5) 



Using (3 

night 



1.2, one 

day , a very significant variation of the 



10, V T ~ 0.15 and T d ay/Tni g ht 

condensation pressure. This implies that air flowing on 
constant pressure levels around the planet would lead to 
a rapid depletion of any condensing species on the day 
side, compared to what would be predicted from chemical 
equilibrium calculations. This can potentially also remove 
important absorbing gases from the day side, as in the case 
of TiO, which can be removed by CaTiC>3 condensation, 



or Na, removed by Na2S condensation (Lodders 199E). Of 
course, most of the variation depends on the exponential 
factor /?/(/? Vt — 1), which is infinite in the limit when 
the atmospheric temperature profile and the condensation 
profile are parallel to each other. 

In the discussion, we implicitly assumed /3Vt — 1 > 0; 
however, when the atmosphere is close to an isotherm, this 
factor can become negative. In this case the day/night 
effect is even more severe, as the condensing species is 
entirely removed from this quasi-isothermal region. 
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